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ABSTRACT 

Context, hya emission is commonly used as star formation tracer in cosmological studies. Nevertheless, resonant scattering strongly 
affects the resulting luminosity, leading to variable and unpredictable escape fractions in different objects. 

Aims. To understand how the Lyo- escape fraction depends on the properties of the star-forming regions, we need high spatial resolution 
multiwavelength studies of nearby Lyo- emitters, like Haro 2. 

Methods. We study the Lyo- emission of Haro 2 in connection with the properties of the young stellar population, the characteristics 
of the interstellar medium, the distribution and intensity of the Balmer emission lines and the properties of the X-ray emission. We 
have used HSTFOC UV, WFPC-2 optical and NICMOS NIR images, as well as STIS spectral images along the major and minor axes 
of Haro 2 to characterize the Lyo- emission and analyze the properties of the stellar population. WFPC-2 Ho- image and ground-based 
spectroscopy allow us to study the Balmer emission lines. Finally, Chandra/ACIS X-ray images provide resolved distribution of the 
X-ray emission at various energy bands. The observational data are compared with the predictions from evolutionary synthesis models 
to constrain the properties of the star formation episode. 

Results. The UV, Ho- and far infrared luminosities of the Haro 2 nuclear starburst are well reproduced assuming a young stellar 
population with ages ~3.5 - 5.0 Myr, affected by differential intestellar extinctions. Whereas the diffuse soft X-ray emission observed 
is attributed to gas heated by the mechanical energy released by the starburst, an unresolved component (likely an UltraLuminous X- 
ray source) accounts for the hard X-ray emission. Both compact and diffuse Lyo- emission components are observed along both axes. 
Lyor is spatially decoupled from Balmer lines emission, Balmer decrement and UV continuum. However, the diffuse Lyo- component is 
spatially correlated with the diffuse soft X-ray emission. Moreover, unlike the compact Lyo- emission, diffuse Lyo- shows luminosities 
larger than predicted from Ha, assuming case B recombination and considering the dust extinction as derived from Ho-/H/3. 
Conclusions. The presence of outflows and dust abundance strongly affects the Lyo- emission closely associated to the massive stellar 
clusters, leading to even a range of escape fractions at different locations within the same starburst. On the other hand, we propose 
that the diffuse Lyo- emission originates in gas ionized by the hot plasma responsible for the soft X-ray radiation, as suggested by 
their spatial correlation and by the measured L(Ha)/Lo.4-2.4keV ratios. Calibration of hya as star formation rate tracer should therefore 
include both effects (destruction vs. enhancement) to avoid biases in the study of galaxies at cosmological distances. 

Key words. Galaxies: starburst - Galaxies: star formation - Galaxies: ISM - Ultraviolet: galaxies - Cosmology: observations - 
Galaxies: individual: Haro 2, Mrk 33, Arp 233, UGC 5720 



1. Introduction 

The spectral energy distribution (SED) of galaxies experienc- 
ing intense and relatively short episodes of massive star forma- 
tion (starbursts) is dominated in the UV by the continuum of 
the young, most massive stars. The ionizing radiation they pro- 
duce interacts with the surrounding gas to produce nebular emis- 
sion lines through recombination processes. The mechanical en- 
ergy released during the evolution of these stars, in the form 
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of stellar winds and supernova explosions, heats the medium, 
producing X-ray emission and collisionally excited emission 
lines. While Ha is the most prominent hydrogen nebular line 
in the optical, simple nebular models predict Lya to be ~9 times 
more intense than Ha when considering typi cal nebular parame - 
ters for temperature and electronic density dOsterbrockl W 9 8 9 ) . 
Nevertheless, observations over the last decades have shown 
that the value of the L(Lya)/L(Ha) may indeed be very dif- 
ferent than expected, and that the presence, intensity, spectral 
shape and spatial profile of Lya line do depend on many fac- 
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tors (s e elMas-Hesse et alJ(l20Q3h . lAtek et alJ (l2009h . lHaves et alJ 
(2010), lOstlin et alJ (l2009h" and references therein). Lya is a 
resonant line and this causes Lya photons to be absorbed and 
reemited by neutral hydrogen atoms within the gas many times 
before they can escape. The increase in the path of the pho- 
tons at Lya wavelength due to this resonant effect implies that 
very little dust can completely destroy the photons around 1216 
A, yielding a d amped Lya absorption dChen & Neufeldlll994t 
lAtek et alJl2009h. However, as explained by iMas-Hesse et alJ 
(120031) and IVerhamme. Schaerer & Masellil (120061) . when the 
mechanical energy injected by the starburst into the interstellar 
medium (ISM) accelerates the surrounding neutral gas, form- 
ing an expanding superbubble, only the blue wing of the emis- 
sion line is scattered. As a result, Lya emission can escape 
from the region, but will show a char a cterist ic P Cyg emission- 
absorption profile. Mas -Hesse et alJ (120031) showed different 
cases in a sample of blue compact galaxies in the Local Universe, 
ranging from damped absorption to emission with and with- 
out P Cyg profiles, depending on the kinematical and ionization 
properties of the medium surrounding the star-forming regions. 
IVerhamme. Schaerer & Masellil (120061) computed the different 
profiles expected applying radiation transfer models to differ- 
ent scenarios, including the presence/absence of accelerated gas, 
du st, energy shift of the photo ns by the scattering process,... 

iTenorio-Ta gle et al J (11999) proposed a sequence of Lya pro- 
files following the different evolutionary phases of the starburst, 
assuming the acceleration of the neutral gas by the ongoing star- 
burst and the formation of expanding shells to be the driver of the 
Lya line visibility. Since the acceleration of the surrounding in- 
terstellar medium is responsible for both t he generation of X-ray 
emiss ion, by the process of gas heating dStrickland & Stevens! 
Il999h . and for the escape of the Lya emission line, by opening 
kinematical 'holes' in the neutral gas, we expect some intrinsic 
correlation between both phenomena. In order to study this pos- 
sible correlation we embarked in a project to study the X-ray 
properties of Lya emitters in the Local Universe, starting with 
Haro 2. 

Haro 2, also known as Mrk 33, Arp 2 33 and UGC 5 720, 
is a metal-rich blue compact dwarf galaxy (iDavidgdl 19891) . ex- 
perie ncing an intense star-form ing episode with an age of ~5 
Mvr dMas-Hesse & Kunthll999h . Star formation is located in the 
nucleus of the galaxy where it can be resolved into individual 
knots (Fig. [TJ. Chandar, Leit herer & Tremontil J20 04) combined 
the HST/STIS spectra of six different clusters with in the star- 
bursting nucleus, whereas M endez & Estebanl ((2000) could only 
resolve three knots in an Ha image obtained with the HIRAC 
camera of the Nordic Optical Telesco pe. On the other hand, 
Summers. Stevens & Strickland! ((2001) detected two knots in a 
B-band image from the Jacobus Kapteyn Telescope. Besides 
the youn g population produced in t he cu rrent star-formation 
episode Fanelli, O'Connell & Thuanl d!988l) argued for the ex- 
istence of an older stellar population given the spectral features 
of A-type stars found by t he International Ultraviolet Explore r 
(IUE). Based on the data bvlFanelli, O'Connell & Thuanl d 19881) . 



Sum mers. Stevens & Strickland (2001) decomposed the star for- 
mation history of Haro 2 into a current starburst of ~5 Myr, 
a previous burst of ~20 Myr and an earlier episode having 
ocurred ~500 Myr ago . Desp ite being a metal rich, rather dusty 
galaxy, Leaue ux et al.1 {l995) found that Haro 2 is an intense 
Lya-emitter. The line shows a clear P Cyg profile apparently 
originated by a neutral superbubbl e expanding at ~200 km s , 
energized by the central starburst. IMas-Hesse et al.1 (120031) ob- 
tained a high-resolution spectral image with H5T/STIS in the 
Lya range along the minor axis, showing the P Cyg profile and 



the spatial decoupling of the continuu m and the Lya emission. 
ISummers. Stevens & Strickland! ([2001) performed a study of the 
X-ray radiation of Haro 2 based on ROSAT/HR1 data and proved 
that the mechanical energy released by the nuclear starburst was 
indeed powering the diffuse, soft X-ray emission. These authors 
quantified in 2% the fraction of the mechanical energy injected 
by the burst into the medium ending up being emitted as X-rays. 
iLegrand et al.l d 19971) identified the large Ha shell at the edges of 
the expanding superbubble, decoupled from the rotation velocity 
of the galaxy. 

The bolometric luminosity of Haro 2 is well represented by 
its infrared luminosity Lfir= 1.4 x 10 43 erg s _1 , which overlaps 
with the Lboi low-end in the luminosity fu nction of Lyman Alpha 
Emitting Galaxies (LAEs) at z = 3.1 dGronwall et al.ll2007l) . 
Haro 2 also shows an UV continuum luminosity similar to that 
of the weakest Lyman Break Galax ies (LBGs) in the redshift 
range z = 2.7-5 dTapken et al.l2007l) . This fact, together with its 
high Lya-intensity and P Cyg spectral profile, makes of Haro 2 
a good local prototype of high-z Lya-emitters whose star forma- 
tion activity is usually derived from both Lya luminosity, UV 
continuum and/or X-ray emission. Understanding the processes 
determining these emissions in a closeby galaxy which can be 
studied in high detail is of paramount importance to quantify the 
star formation activity in distant galaxies. The basic properties 
of Haro 2 are summarized in Table Q] 

In this paper we have analyzed observations obtained by HST 
with STIS, NICMOS and WFPC-2, and by Chandra in X-rays, 
as well as data obtained from ground-based observatories, to 
fully characterize the properties of the star formation episodes 
taking place in the nucleus of Haro 2 and of its associated X-ray 
and Lya emissions. In Sect. [2] we describe the HST, Chandra 
and ground-based observations. In Sect. [3] we present the results 
obtained from the analysis of the observational data, which we 
discuss later in Sect. [4] Sect. [5] contains the conclusions of the 
study. 

Throughout this paper we will consider 20.5 Mpc as the 
distance to Haro 2, which yields a projected scale of ~100 pc 
arcsec -1 . We want nevertheless to remark that the distance to 
Haro 2 could be larger, when compared to IZwl8, whose es- 
timated distance has been increased by a factor of almost 2 to 
19.5 Mpc by the analysis of Colour Magnitude Diagrams and 
Cepheid variables. The absolute values concerning luminosities, 
masses and so on should be taken with a word of caution, al- 
though they do not affect our conclusions. 



2. Observational data 

2.1. HST observations: ultraviolet, optical and near infrared 

We have used observations of Haro 2 obtained with STIS, 
WFPC-2 and NICMOS, in the UV, optical and near infrared 
ranges, obtained from the HST Archive at STScI. In this section 
we describe the data and the processing done. 

2.1.1. HST images 

The HST imaging data sets used are listed in Table|2] There is a 
single UV image of Haro 2, obtained with FOC in 1993, before 
COSTAR was installed. The observation was performed with the 
F220W filter. The far-UV image covers -21" x 23" with plate 
scale of ~0.02" pixel , and is centered on A = 2328 A. The 
FOC image is shown in Fig. [2] together with the STIS slits direc- 
tions perpendicular to the major axis (direction Dl) and parallel 
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Fig. 1. Left: image of the whole galaxy Haro 2 from NOT/ ALFOSC in the V-band. Right: image of the UV-continuum of the nucleus 
of Haro 2 obtained with HSTfWFPC-2 with filter F336W where knots SE and NW are labelled. 

Table 1. Haro 2 coordinates, distance, scale, Galactic HI, color excess towards the source, and oxygen abundance. 



R. A. 


Decl. 


Redshift" 


Distance a 


Scale 3 


yv(Hl) b 


E(B-V)aa 


12 + log(0/H) c 


(J2000.0) 


(J2000.0) 




(Mpc) 


(pc arcsec -1 ) 


(cm" 2 ) 






10 32 31.9 


+54 24 03.7 


0.004769 


20.5 


100 


6.3 x 10 19 


0.012 


8.5 



a From NED (the NASA/IPAC Ext ragalactic Database). 
b Value from lLeaueux et all d 19951) . 
c Value from Davidgel dl989l) . 



to it (direction D2). Since UV emission is dominated in star- 
forming regions by massive stars, with this FOC image we can 
study the UV-morphology of the starburts located in the nucleus 
of the galaxy, as well as calculate the UV-flux at A — 2328 A 
(Fuy, hereinafter) emitted by each knot and by the total source 
as a whole. To derive the Fyv f ux of each knot we defined two 
similar and non-overlapping boxes, as indicated in Fig. [2] The 
boundaries of these boxes are rather arbitrary, since we can not 
discriminate beween both components, but the flux inside them 
should be dominated by the bright stellar clusters within knots 
NW and SE. The sizes of the boxes are 7.4" x 5.5" (knot NW) 
and 7.4" x 4.8" (knot SE). We have checked that the brightest 
knots account indeed for 70 - 80% of the total UV flux measured 
on each box, dominating so their integrated spectral properties. 
As we will discuss later, while the properties of the starburst will 
be therefore dominated by the massive stars in the bright knots, 
we have used the total flux within the boxes to derive the mass 
normalization when applying evolutionary synthesis models. 

Background on the FOC image was estimated from the most 
separated regions away from the source, averaging on boxes of 
3 x 3, 5 x 5 and 10 x 10 pixel 2 . The UV-background was esti- 
mated to be ~0.75 counts pixel" 1 . In order to work with an image 
free from background contamination, we subtracted this average 
background from all pixels of the image. After that, the intensity 
was multiplied by PHOTFLAM and divided by the total exposure 



time EXPTIME in order to work with fluxes in units of erg s 
cm" 2 A" 1 . 

We studied the effect of FOC point spread function (PSF) 
since the image was obtained prior to COSTAR deployment. We 
downloaded a pie-COSTAR PSF for an F/96 F220W observation 
performed five months earlier from the instrument websitf] and 
extracted the enclosed point source flux within the boundaries 
of the boxes we have considered. It was found that for the sizes 
and locations of the regions considered and assuming the knots 
are point-like sources, a ~85% of the total radiation is included 
within the integration area, hence leading to an underestimation 
in the total Fuv of less than ~15%. 

Haro 2 has been observed with H5T/WFPC2 in the UV- 
optical range with broad-band filters F336W, F439W, F555W 
and F814W, providing images of the continuum emission, and 
with the narrow filter F658N, which is dominated by Ha. In or- 
der to isolate the Ha-emitting regions, continuum was subtracted 
from the F658N image using the F814W one, assuming it rep- 
resents well the continuum at Ha. After downloading the data 
from the HST archive, drizzled files were combined and aligned 
for filters F814W and F658N. Then, the resulting images were 
multiplied by PHOTFLAM in order to convert counts into flux val- 
ues. Finally, the F814W image was multiplied by an arbitrary 



1 http://www.stsci.edu/hst/foc/calibration/f96_novl992.html 
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Fig. 2. L eft: Ultraviolet HST/FO C image of Haro 2, together with the locations of the HST/STIS s lits Dl (gratings G140M and 
G430L) dMas-Hesse et alJ2003l) . and D2 (G140L and G230L) dChandar. Leitherer & Tremontill2004h . and the boxes used to extract 
the UV flux from each knot. The width of the slits is represented at real scale. North is up and East is left. Right: Ha image obtained 
with HST/WFPC2, with the UV continuum contours superimposed. 



constant and then subtracted from the F658N image, using field 
stars as reference. 

Finally, we downloaded from HST archive an image of 
Haro 2 obtained by NICMOS at 1.60, ^m with broad-band filter 
F160W. Its size is 19.2" x 19.2", and the scale of the pixel is 
0.075". No background subtraction on the science data file was 
necessary since NICMOS reduction pipeline already performs 
this removal through task cnib. The NICMOS F160W, WFPC- 
2 F336W and WFPC-2 Ha images have been colour-coded and 
combined in Fig. [3] 

The astrometry included in the World Coordinate System 
(WCS) keywords in the headers of the different images was not 
always consistent, with relative offsets of up to ~1.2" between 
the NICMOS and FOC images, for example. We have taken the 
HST/NICMOS header information as reference, which matched 
well the Chandra one, and have corrected the astrometry from 
other images when needed. 

2.1.2. HST spectroscopy 

M as-Hesse et al.l ((2003) studied spectral images obtained with 
the G140M and G430L grating on STIS, with a 52" x 0.5" long 
slit placed orthogonal to the major axis of the nucleus. This high 
resolution spectral images covered the range 1200 - 1250 A, 
including the Lya line, and 2900 - 5700 A, comprising the oxy- 
gen lines 0[II]3727 and 0[III]4959, 5007, and H/3. Although the 
integration time on this optical spectral image was very short, 
it provides the spatial profiles of these lines with high angu- 
lar resolution. Hereinaft er, we will refer to the direction of this 
slit as Dl. Furthermore. IChandar. Leitherer & Tremontil 12004) 
analyzed two additional spectral images using the G140L and 
G230L gratings with the 52" x 0.2" long slit oriented along 
the major axis of the starbursting region. These spectral im- 
ages cover the spectral ranges 1 150 - 1730 A, showing the Lya 
emission and the stellar absorption lines Si IV 1394, 1403 and 



C IV 1548, 1551, and 1570 - 3180 A. This direction will be re- 
ferred to as D2. Information on both UV-optical observations is 
listed in Table[3] and in Fig.|2]both slit directions can be observed 
superimposed on the FOC image. Both sli t directions are or- 
thogonal, and whereas livIas-Hesse et al](l2003l) encloses only the 
SE sta r-forming knot, slit from Chandar, Leitherer & Tremonti 
d2004l) contains the bulk emission from both SE and NW knots. 
The G140M Dl and G140L D2 spectral images constitute the 
direct observation of the Lya morphology and spectroscopy of 
Haro 2 which we analyze in this work. 

Extraction, background removal and anal ysis of the S TIS 
Dl da ta were already discussed in detail by iMas-Hesse et al.l 
(I2003h . Data were reduced with STSDAS procedures and flux 
from 100 rows above and below the stellar continuum were av- 
eraged to estimate the geocoronal Lya emission and the detector 
background, being the latter negligible. Slit correction was per- 
formed by multipliying measured counts by the STIS parameter 
diff2pt. Observational data for filters G140L and G230L in 
D2 were downloaded from the HST archive already calibrated. 
Again, intensity values measured were multiplied by diff2pt 
in order to correct for slit losses. After that, the first step was 
to select the continuum regions, deriving continuum extensions 
within the slit of ~60 pc and ~ 150 pc for both knots SE and NW 
respectively (see Table [6]). The integration of the contributions 
from each of the regions of the spectral images, and the merging 
of the observations performed with each grating were done with 
IRAF-STSDAS procedures. 

We noticed that keywords RA_APER and DEC_APER in the 
headers of the G140M and G430L HST/STIS file did not corre- 
spond to the actual position of the slit, centered on knot SE, but 
rather referred to an intermediate position closer to knot NW. 
From the Ha image obtained with HSTfWFPC-2 data we ex- 
tracted two different spatial profiles, one assuming the position 
of the slit indicated by the HST/STIS observation file and which 
is centered close to knot NW, and another profile for a slit placed 
on knot SE. The angle of the slit was fixed to the value indicated 
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Fig. 3. Multiwavelength images of Haro 2: NICMOS F160W near-infrared continuum (red), WFPC-2 F336W UV continuum (blue) 
and Ha emission (green). Images are 10 arcsec wide. The black cross marks the center of the bluest stellar cluster in knot SE. Note 
the red stellar population to the East of knot SE, and disperse within knot NW. North is up and East is left. 

Table 2. Log of HST/FOC, HST/NICMOS, HST/WFPC2 and Chandra observations of Haro 2 used in this work. 



Instrument 


Observation date 


Filters 


Integration time 

(s) 


/!o 
(A) 


Plate scale 
(arcsec pixeP 1 ) 


HST/FOC* 


1993 Apr 18 


F220W 


1197 


2328 


0.02 


//S77NICMOS 


2006 Sep 28 


F160W 


912 


1.6 x 10 4 


0.075 


tfST/WFPC2 b 


2009 Mar 9 


F336W 


3300 


3332 


0.10 


#S77WFPC2 b 


2009 Mar 9 


F658N 


1800 


6591 


0.10 


tfST/WFPC2 b 


2009 Mar 9 


F555W 


600 


5304 


0.10 


#S77WFPC2 b 


2009 Mar 9 


F814W 


600 


8337 


0.10 


Chandra/AClS c 


2008 Feb 2 




19194 




0.49 



Observation performed before COSTAR deployment. 



in the file header. Given the same nature of Balmer lines Ha and 
HjS, we compared the extracted Ha profiles with the H/3 profile 
from HST/STIS G430L observation, and they only agree when 
the slit is assumed to be centered on knot SE. We did the same 
comparison between profiles extracted from the HST/WFPC-2 
F555W image and the profile of the composite H/?+0[III]4959, 
5007 +continuum obtained from the HST/STIS G430L observa- 
tion and the match is only possible when the slit is placed on 
knot SE. Therefore, we concluded that HST/STIS G140M and 
G430L observations of Haro 2 are centered on knot SE and in 
Fig.[2]this slit is located accordingly. 



2.2. Chandra: X-rays 

Chandra time was allocated under proposal 09610593 (PI: J. 
Miguel Mas-Hesse) to observe Haro 2 with ACIS. Observation 
9519 with time duration of ~20 ks (see Table 13 was performed 
on February 2, 2008. Chandra/ 'ACIS data were reduced with 
CI AO 4.1 following the standard reduction procedures available 
on the Chandra-CIAO website, reprocessing level 1 and 2 event 
files. The Ancillary Response File (ARF) and the Redistribution 
Matrix File (RMF) were also calculated. In order to obtain the 
X-ray spectrum of Haro 2 source counts were extracted from 
a circular region of radius ~9" centered on the bulge emisison 
since it incloses most of the X-ray radiation of the source. This 
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Table 3. Log of HST/STIS observations of Haro 2 used in this work. 



Observation date 


Grating 


Integration time 


Position angle 


Slit direction 


Slit size 


Wavelength interval 


Plate scale 


Spectral dispersion 






(s) 


(deg) 






(A) 


(arcsec pixer 1 ) 


(A pixer 1 ) 


2000 Dec 1 


G140M 


7900 


-121.35 


Dl 


52" X 0.5" 


1200- 1250 


0.029 


0.053 


2000 Dec 1 


G430L 


300 


-121.35 


Dl 


52" X 0.5" 


2900 - 5700 


0.050 


2.746 


2003 Feb 26 


G140L 


1446 


146.85 


D2 


52" X 0.2" 


1150- 1730 


0.025 


0.584 


2003 Feb 26 


G230L 


600 


146.85 


D2 


52" X 0.2" 


1570-3180 


0.025 


1.548 



Table 4. Log of WHT/ISIS observations of Haro 2 used in this work. 



Observation date 


Arm 


Grating 


Integration time 


Position angle 


Slit direction 


Slit size 


Wavelength interval 


Plate scale 


Spectral dispersion 








(s) 


(deg) 






(A) 


(arcsec pixer 1 ) 


(A pixel" 1 ) 


2006 Dec 25 


Blue arm 


R600B 


3600 


-121.35 


Dl 


3.3' X 1.0" 


3460 - 5280 


0.20 


0.45 


2006 Dec 25 


Red arm 


R600R 


3600 


-121.35 


Dl 


3.3' X 1.0" 


5690 - 7745 


0.22 


0.49 


2002 Feb 9 


Blue arm 


R300B 


3600 


146.85 


D2 


4.0' X 1.1" 


3480 - 7020 


0.20 


0.86 




region is shown in Fig. [4] labeled as NUCLEUS. Background 
was extracted from a composite of three circular regions with 
radius ~46", 56", 58" from the area surrounding the galaxy and 
not containing any other bright source. A net total count rate of 
~1.4 x 10~ 2 cts s _1 after background subtraction was measured 
for the whole energy range in NUCLEUS. Sp ectral analysis was 
performed with XSPEC v. 12.7 (lArnaudlll996h . grouping source 
counts to have at least 15 counts per bin in order to be able to 
a PPly X 2 statistics in the fitting analysis. The final binned X-ray 
spectrum of Haro 2 is shown in Fig. [8] 

Chandra/ AOS X-ray image is displayed in Fig. [4] colour- 
coded as a function of the energy, where red: 0.2-1 .5 keV, green: 
1.5 — 2.5 keV and blue: 2.5 - 8.0 keV. Color-coded X-ray image 
is also shown in Fig.|5]with higher detail, together with the com- 
posite multiwavelength image at the same scale for reference. 
We have also superimposed the NICMOS near-infrared contours 
to identify the positions of the star-forming knots. Haro 2 shows 
a diffuse soft X-ray emission extended mostly to the N and W 
of the starbursts, over a region with radius ~600 pc. Also, three 
hard, point-like sources are observed in Fig. [4] and labeled as 
XI, X2 and X3. Background-corrected count rates of these hard 
sources for several energy bands are shown in Table |5j together 
with the values for NUCLEUS, which includes XI and X2. XI 
is located at the position of knot SE, and is associated to its 
massive stellar cluster, as we will discuss later. X2 is a source 
5 times weaker in the 1.5 - 8.0 keV band than XI, located to 
the Southeast of knot SE. X2 does not have any correlation to 
any continuum or Her emitting region, since, as Fig.|5]shows, the 
closer UV-optical bright knot is ~1.5" far. An extra hard object 
X3 detached from the starburst region is found ~10" Northwest 
from the nucleus, with a net flux of ~ 1 .5 x 10~ 3 cts s _1 . The spec- 
tral image obtained with NOT/ISIS in direction D2 (see Sect. 12. 3t 
indicates that the ratio 0[II]3727/H/3 in the region close to X3 is 
two orders of magnitude higher than in the nuclear region. This 
fact points towards X3 being a source within Haro 2 rather than 
a background source. However, given its angular distance from 
the bulk of the X-ray emission and its low flux, we did not an- 
alyze it in this work. Therefore, XI is the dominant component 
of the hard emission in the NUCLEUS. 



2.3. Ground-based observations 

An image of Haro 2 obtained with instrument ALFOSC on 
the Nordic Optical Telescope with filter V_Bes 530_80 is dis- 



Fig. 4. Color-coded Chandra X-ray image of Haro 2: soft (0.2 - 
1.5 keV, red), medium (1.5 - 2.5 keV, green) and hard X-rays 
(2.5 - 8.0 keV, blue). The contours of the NICMOS near-infrared 
continuum and the position of the STIS slits have been super- 
imposed. Main X-ray components are labeled. Source region 
from which X-ray spectrum was extracted is confined within the 
dashed circle. Regions within the yellow ellipses show diffuse 
Lya emission. 



played in Fig. Q] (courtesy of Luz Marina Cairds and Ricardo 
Amorin). Also, two long-slit observations with ISIS on the 
William Herschel Telescope were downloaded from the archive 
for which information is included in Table [4] ISIS spectral im- 
ages of the source were acquired with slit aligned along direc- 
tions Dl (red and blue arm in December 25, 2006) and D2 (blue 
arm in February 9, 2002). Although the observation along D2 
was performed only with the blue arm of the instrument, a sim- 
ilar spectral range than the one achieved for direction Dl with 
both arms was obtained using an observational mode with lower 
spectral dispersion. As a result, spatial configuration of stellar 
continuum and main nebular emission lines was obtained along 
both the major and minor axes of Haro 2. Specifically, we de- 
rived from these images the spatial profiles of the Hy8 and Ha 
emission lines. 
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2 arcsec 




Fig. 5. Left: Multiwavelength image of Haro 2, as in Fig. [3] Right: the color-coded Chandra/ACIS X-ray image (colors as in Fig. |4| 
at the same scale, with the NICMOS contours superimposed for reference. North is up and East is left. 

Table 5. Count rates of X-ray objects corrected for background contamination. 



Object 


Total 

(s- 1 ) 


Soft(0.2- 1.5 keV) 

(s- 1 ) 


Medium(1.5 - 2.5 keV) + Hard(2.5 - 8.0 keV) 

(s" 1 ) 


NUCLEUS 


2.o x i(r 2 


1.2 x itr 2 


3.2 x 10~ 3 


XI 


6.8 x 1(T 3 


4.2 x itr 3 


2.1 x 10~ 3 


X2 


8.9 x 1(T 4 


3.6 x itr 4 


4.2 x 10~ 4 


X3 


1.5 x 10" 3 


4.2 x itr 4 


6.3 x 10~ 4 



3. Analysis and results 

3.1. HST morphology and photometry 

The UV image of Haro 2 nucleus with the identification of its 
brightest star-forming knots is shown in Fig. [TJ together with 
an image of the whole galaxy in the V-band. The elliptically- 
shaped nucleus of Haro 2 ho sts various massive cluster s, lo- 
cated on knots SE and NW. iMendez & Estebanl d2000l) pro- 
posed to differentiate 2 smalle r knots within NW, as indeed 
sugges ted by Fig. [2] Moreover, Chandar, Leith erer & Trem onti 
(2004) argued they could identify 6 faint clusters within the 
galaxy, from the inspection of the UV spatial profile along slit 
D2 (see Fig. |7). However, for the purposes of this work we will 
consider only knots NW and SE as individual star-forming re- 
gions. These knots are aligned along the major axis, being sep- 
arated by -75 - 100 pc. While knot SE is not resolved, knot 
NW is more extended and shows a more complex morphology. 
UV continuum extends few hundreds of parsec away from the 
knots, so that the whole nuclear star-forming region has a size of 
-750 x 1000 pc 2 . 

The UV-flux values for each knot and for the whole star- 
burst region were extracted from the FOC image considering 
boxes shown in Fig.|2]with sizes 7.4" x 5.5" for knot NW and 
7.4" x 4.8" for knot SE. As explained above, due to the proxim- 
ity between knots NW and SE the origin of the counts could not 
be discerned unambiguously, hence the total flux of the whole 



source should be taken as a more reliable measure. Fuv values 
for each knot together with their sizes are included in Table [6] 
The total observed UV emission of the whole galaxy nucleus 
detected by FOC is Fw= 2.7 x 10~ 14 erg s" 1 cm" 2 A" 1 which 
agrees with the value from the spectrum taken by IUE with a 
total aperture of 10" x 20", which should enclose the whole UV- 
emitting region. As it will be explained later, masses of the bursts 
were calculated through the UV flux of each of the knots using 
evolutionary population synthesis models. 

Comparison of NICMOS image of IR continuum at 1.60 
fim with the WFPC-2 image of the UV continuum at 3300 A 
(see Fig. [3]l shows that some regions which are very prominent 
in the IR are not visible in the UV. This results especially evi- 
dent on knot SE, as well as in some more diffuse areas of knot 
NW. Since both knots are dominated by young, massive stars, 
the non-detection in the UV of the stars detected in the IR is 
probably due to a patchy distribution of dust, which hides part 
of the young population, or to evolutionary effects on different 
populations. 

The Ho- image obtained with HST/WFPC2 is shown in Fig. [2] 
(right), with the UV continuum contours superimposed. The Ha 
emission appears associated and essentially co-spatial to both 
stellar clusters SE and NW. 
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Table 6. Physical properties of knots NW and SE in Haro 2, assuming a single stellar population at each knot, with mass normal- 
ization according to their integrated UV luminosity (see text for details). 



Knot 


Age 


Mass 


E(B-Vf 






Stellar extension c 




(Myr) 


(M Q ) 




(erg s _I cnT 2 A -1 ) 


(erg s _I A -1 ) 


(pc) 


Knot SE 


~4 


7.0 x 10 s 


0.050 


1.1 x 10~ 14 


5.5 x 10 38 


-60 


KnotNW 


~5 


1.3 x 10 6 


0.020 


1.7 x 10 14 


8.5 x 10 38 


-150 



a Internal color excess assuming Pre vot et al.l (Q984) extinction law (SMC). 

b Luminosity value not corrected for Galactic or intrinsic extinctions. 

c Extension of the stellar continuum along direction D2 in the STIS spectral image. 



Table 7. Predicted values of the ages, masses, stellar and nebular extinctions and emissions of Populations 1, 2 and the composite 
of both, as well as the EW(HB) produced. Observational values are also displayed. 



Population Age M E(B-V), E(B-V) neb L uv " L(Hor) a L FIR EWQiB) 





(Myr) 


(M ) 






(erg s 1 A ') 


(ergs ') 


(ergs ') 


(A) 


1 


4.5 


1.8 x 10" 


0.035 


0.24 


1.3 x 10 39 


2.5 x 10 40 


3.4 x 10 42 




2 


3.5 


2.0 x 10 6 


0.5 


0.24 


7.1 x 10 37 


5.3 x 10 40 


1.6 x 10 43 




1+2 




3.8 x 10 s 




0.24 


1.4 x 10 39 


7.8 x 10 40 


2.0 x 10 43 


58 


Observed 








0.24 b 


1.4 x 10 39 


7.8 x 10 40 b 


2.0 x 10 43 


44 - 80" 



a Luminosit y value reddened by Galactic an d intrinsic extinctions. 
b Value from lMoustakas & Kennicuttl d2006h . 



3.2. UV continuum 

High spe ctral resolution UV s pectra are available in Starburst99 
models (Leith erer et al.1 (Il999l) . hereinafter SB99 models) in the 
range 1 200 - 1 600 A for the LMC/SMC library and 1 200 - 1 800 
A for the Milky Way library, with a spectral step of 0.75 A. 
Normalized spectra for both knots in the wavelength range 
1150 - 1700 A were extracted from STIS spectral image ob- 
tained in the direction D2 with filter G140L. The spatial exten- 
sion of the stellar continuum regions considered was 0.6" and 
1.5" for knots SE and NW, respectively. Given the sentivity of 
the shape and intensity of the absorption lines Si IV 1394, 1403 
and C IV 1548, 1551 to the age of the burst, normalized UV spec- 
tra from both knots were fitted individually with SB99 UV mod- 
els in order to calculate their evolutionary state. SB99 models 
were computed for an initial mass function <f>(m)~m~ 2 - 35 (i.e. a 
Salpeter IMF) and mass limits of 2 - 120 M , assuming instan- 
taneous star formation and Z = 0.008 for the evolutionary tracks 
and high-resolution spectra. This metal licity is the clo sest avail- 
able to the observational one Z~ Z /2 (lDavidgelll989h . UV line 
spectra from the models were either taken from the Milky Way 
or the LMC/SMC libraries. We found that no acceptable fitting 
of the normalized spectrum of either knot was possible with the 
LMC/SMC library whereas a good agreement was reached when 
the Milky Way library was used. We obtained ages of ~4 Myr 
and ~5 Myr for knots SE and NW. Uncertainties in the proce- 
dure may yield a precision in the values of the ages obtained 
of approximately +0.5 Myr. The final fittings of the normalized 
spectra are shown in Fig. [6] 

Reddening in Haro 2 was calculated through the fitting of 
the UV SED of each knot with the SB99 models, given the 
assumptions for IMF and metallicity already explained above. 
SEDs of the models at the ages 4 Myr (for knot SE) and 



5 Myr (for knot NW) were reddened for a certain value of 
the color excess E (B-V) and assuming SMC extinction law 
dPrevot et al.ll 1984 and Ry = 3.1, then they were redshifted and 
finally extinguished for Galactic extinction for th e color excess 
E(B-V )- 0.012 assuming Cardelli extinction law (Cardelli et al. 
119891) . Redshift was obtained from the Gaussian fit of the HB 
line included in the spectrum of the G430L observation, which 
turned out to be z = 0.00461, that is within ~4% the value re- 
ported in NED (NASA/IPAC Extragalactic Database). Iterative 
method for calculating E(B-V) was repeated until a good fit 
was achieved for each knot, obtaining E(B-V)k not s,E~0-050 and 
£ , (B-V)knotNW~0.020. Observed spectra are shown in Fig. [6) to- 
gether with the fitting SEDs of the models reddened. Finally, 
the masses of both knots were calculated dividing the observed 
UV luminosity Luv of each knot measured by HST/FOC by the 
reddened, redshifted, predicted value from the SB99 models, 
which is computed per unit mass, obtaining Mi cnot sE~7 x 10 5 
M Q and Mk no tNW~l-3 x 10 6 M Q . These values are obtained 
when considering IMF mass limits of 2 - 120 M . Were other 
boundaries of the mass range be consider ed, a proper correc- 
tion should be performed to the se values dWilkins et al.l f2008; 
lOti-Floranes & Mas- Hesse1l2010l) . Given the difficulties to dis- 
cern the actual flux of each knot from the total one, the values 
of the masses calculated should be taken as a rough estimate. 
As explained in Sect.2.1.1, this procedure assumes that the stars 
contributing to the diffuse UV continuum within each box, as in- 
dicated in Fig. [2] share the same age as those in the central knots 
dominating the UV spectral features. While this should be con- 
sidered just as a first order approximation, this strategy provides 
a more realistic global mass normalization than just considering 
the flux within the unresolved knots. A correct normalization is 
required to properly compare the predictions with globally inte- 
grated magnitudes, such as the X-ray or far infrared luminosities. 
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Fig. 6. Left: observed UV stellar continuum of knot SE (bottom) and knot NW (top). Right: normalized spectrum of UV stellar 
continuum of knots SE and NW (solid lines) with absorption lines Si IV1394, 1403 and C IV1548, 1551 fitted with SB99 models 
for ages ~4 Myr and -5 Myr (dashed lines), respectively (see text for details). 



Physical properties calculated through evolutionary popula- 
tion synthesis models for knots SE and NW are shown in Table 
[6] Given the small volume of knot SE, which we can not resolve, 
and its mass, Mk not sE > 1 x 10 5 M H , this knot should be consid- 
ered a super stellar cluster dAdamo et alJl20l6T) . 



We want to stress that the results of the fit discussed above 
are based on the UV continuum and spectral features as ob- 
served with STIS. Lfir and L(Ha) emissions provide us with 
a consistency check to verify the validity of the fitting. Once the 
parameters of the models were calculated, we compared their 
predictions of Lfir and L(Ha) with the observational values. 
O bserved Lfir lumin osity was calculated through the relation 
of lHelouetai1(ll985h given the IRAS fluxes Leo pm = 4.7 Jy and 
Lioo/jm = 5.3 Jy. Since this relation yields the infrared emission 
collected in the wavelength range 42.5 - 122.5 yum, it should be 
corrected in order to obtain th e total emission in the complete in- 
frared range. As was found bv lCalzetti et alj d2000h for a sample 
of eight starburst galaxies, the ratio FIR(l - 1000 ^ m)/FIR(40- 
120_//m)~1.75, hence the value obtained from iHelou et"ai] 
(fT985l) relation was multipli ed by 1 .75, yielding Lfir — 2.0 xlO 43 



erg s~' . On the other hand, Moustakas & K ennicutfl {2006) mea- 
sured the integrated intensities and the equivalent widths of the 
Balmer and forbidden lines for more than 400 galaxies using 
drizzled spectroscopic observations with Boiler & Chivens spec- 
trograph on Bok telescope. Intensities and equivalent widths of 
all lines were corrected for Galactic reddening and stellar ab- 
sorption, and also for [Nil] contamination in the case of Ha, ob- 
taining for Haro 2 L(Ha)= 1.6 x 10~ 12 erg s cm 2 . Since their 
extraction aperture was of 50" and the slit was drifted along 40", 
this value is expected to represent well the total flux associated 
to the burst. The same F (Ha) value for Haro 2 was measured 
bv lGil de Paz et aL[( 2003). Moreover, besides the integrated val- 



ues, 



Moustakas & Kennicutt (2006) reported the values of nu- 



clear regions using an aperture of 2.5" x2. 5". Nebular reddening 
calculate d through the ratios of the inte grated Balmer lines val- 
ues from lMoustakas & Kennicuttl (120061) is L(fl-VJ=0.24, which 
is in perfect agreement with t he value E(B-V)-Q.22 reported 
by iMas-Hesse & Kunthl (1 19991) . This value is -0.2 dex larger 
than the values found by continuum SED fitting for knots SE 
and NW. Higher extinction for nebular gas than for stellar con- 
tinuum has been observed in se veral star-forming re gions (see 
iMaiz-Apellaniz et all ll 19981) and lCalzetti eTall dlOOO)). This dis- 



agreement is usually attributed to the fact that young clusters 
clean its surroundings from dust, leaving relatively clear holes 
through which the UV continuum escapes with low attenuation. 
On the other hand, the more extended ionized gas still contains 
a significant amount of dust. 

Lfir and L(H g) were calculated with the evo l utiona ry syn- 
thesis models by ICervifio. Mas-Hesse & Kunthl d2002l) (here- 
inafter CMHK02 models). Lfir is produced by the emission of 
dust, which is assumed to be in thermal equilibrium. The mod- 
els compute the energy absorbed by dust, comparing the intrin- 
sic UV-optical continuum with the observed one, and assumes 
that all the missing energy is reemitted as Lfir. In addition, the 
models assume that a fraction 1 - /=0.3 of the ionizing pho- 
tons is directly absorbed by dust and does not contribute to the 
ionization. Lfir and L(Ha) were estimated with the CMHK02 
models for the knots SE and NW, assuming the values (mass, 
age, reddening, ...) obtained in the previous analysis and shown 
in Table [6] When the values of Lfir and L(Ha) were compared 
to the observational ones we found that the models severely un- 
derestimated both the Lfir and L(Ha) luminosities by factors 7 
and 3, respectively. 

This disagreement points to the existence of an extremely 
reddened stellar population, whose contribution to the UV con- 
tinuum is negligible in the FOC image and the STIS spectra, 
but which contributes significantly to both the observed FIR and 
Ha emissions. Moreover, the high underestimation of L(Ha) 
points to a young population, with a relatively strong ionizing 
power. Therefore, in order to explain all observables we built 
a simplified model based on the existence of two populations: 
Population 1, which appears as a conspicuous UV-emitter, and 
Population 2, which is highly reddened and is the main con- 
tributor to Lfir and L(Ha). We defined three equations with 
the observables UV, FIR and Ha fluxes being the sum of the 
reddened contributions by each of the populations, where their 
masses Mi and M2 and the stellar extinction of Population 2 
E(B-V)*2 are the unknowns. The contribution by each popu- 
lation was calculated with the CMHK02 models assuming the 
corresponding parameter values. A mean age of Age(l) = 4.5 
Myr and stellar reddening E(B-V)*,\ = 0.035 were assumed for 
Population 1, combining the results for both knots NW and SE 
and mean nebular extinction was assumed to be E(B-V) ne b- 0.24 
dMas-Hesse & Kunfh|[T999HMotistakas & Kennicudl2006h irre- 
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spective of the population. For each age of Population 2 consid- 
ered (Age(2)), a system of three equations with three unkowns 
was obtained and solved. The results show that in order to re- 
produce the observed values, stars from Population 2 must be 
in a less evolved state than those of Population 1, which agrees 
with the former being more reddened having had less time to 
clean the surroundings. Values Age(2) = 3.5 Myr and stellar 
color excess E(B-V),,2 = 0.5 can account for the total UV, Ha 
and Lfir emissions. Assuming lower values for E(B-V)*i im- 
plies a higher contribution by Population 2 to the UV emission, 
which as we have discussed is ruled out by the FOC image and 
the STIS spectrum. Values found for M u M 2 and E(B-V) % i- 2 
are included in Table [7] together with the ages considered, the 
contribution by each population to the magnitudes observed and 
the global EW(HB). Values from Table[6]do differ from those of 
Population 1 of Table [7] since different assumptions were made 
in each case. The values given in Table[6]correspond to the aver- 
age properties of the subsample of massive stars which are less 
affected by interestellar extinction, while a similar fraction of 
young stars seems to be still too obscured to be detected in the 
UV. While this is a simplified scenario, it shows that in reality 
the knots of Haro 2 are probably formed by young massive stars 
with an age within ~3.5-5 Myr, and a continuous sequence of 
interstellar extinction, spanning from regions almost devoid of 
dust to others mostly blocked in the UV by dense dust clouds. 

The underestimation of the Balmer and FIR emissions in 
Haro 2 when assuming a starburst mass derived from the UV 
flux was also found bv lMas-Hesse & Kunthl (1 19991) who, assum- 
ing a single reddened stellar population, obtained lower L(HB) 
and Lfir values than obser ved by factors 2.3 an d 2. Th ese fac- 
tors differ from ours since iMas-Hesse & Kunthl dl999h consid- 
ered a steeper IMF with slope a = -3. On the contrary, they 
obtained an age for the whole nuclear star-forming region of 4.8 
Myr using EW(HB) corrected for both reddening and underly- 
ing stellar absorption, which agr ees with our estimated value- 
Using an updated set of models ICervino. Mas-Hesse & Kunthl 
(2002) found a value of 4. 2 Myr assuming the ver y 
same value of EW(UB) from IMas-Hesse & Kunthl dl999l) . 
IChandar. Leitherer & Tremontil (120041) fitted the UV spectrum 
from G140L observation along direction D2 using SB99 models 
and giving different weights to continuum, interstellar lines and 
stellar wind lines, obtaining the best fit for an age of 5 ± 1 Myr. In 
all these c ases the study was perfo r med o ver the whole galaxy. 
Moreover. iMoustakas & Kennicutil d2006l) obtained EW(BB) = 
44 A and EW(HB) = 80 A for the whole galaxy and the nu- 
clear region, respectively. Both values are already corrected for 
the Galactic and internal extinctions and stellar absorption. Our 
prediction from the models when assuming the existence of 
Populations 1 and 2 of EW(HB)~58 A lies between both values. 

3.3. Lya emission 

Spectral images on the ultraviolet along the minor (Dl direc- 
tion) and major (D2 direction) axes in the Lya wavelength range 
are displayed in Fig. [7] and a relatively strong Lya emission is 
observed in both cases, extending well beyond the stellar clus- 
ters. In the case of the observation of knot SE along Dl direc- 
tion the Lya-emitting region extends largely to the southwest 
directio n whereas it is much w eaker Northeast of the stellar 
cluster. IMas-Hesse et al.l d2003l) showed that the spatial exten- 
sion of this Lya emission along this direction is much larger 
than the continuum-emitting region. Even more, as observed 
in Fig. [7] there seems to be a strong, compact component near 



the stellar continuum, which extends over ~250 pc (Lya;, here- 
inafter, F(Lya) = 2.0 x 10~ 13 erg s _I cirT 2 not corrected for any 
kind of absorption), and another weaker and more diffuse emis- 
sion over ~850 pc extending to the southwest direction (Lyajj, 
F(Lya) = 2.8 x 10~ 14 erg s _1 cnr 2 ). 

In the spectral image along D2 direction displayed in Fig. [7] 
Lya emission seems to be more complex, with regions of total 
absorption and total emission. Fig. [7] shows the spatial profile 
along direction D2 of Lya, being its value positive whenever it 
appears as emission and negative were it pure absorption. The 
UV stellar continuum from HST/STIS is superimposed to the 
Lya-profile and a clear spatial detachment is observed between 
both emissions. Whereas the stellar continuum shows two peaks, 
one rather narrow due to the knot SE (~60 pc) and the other 
one more extended owing to knot NW (~150 pc), Lya-profile 
is found in emission in three different regions: 1) Northwest of 
knot NW, totally detached from the stellar continuum and ex- 
tending over > 600 pc (emission Lya;;;, F(Lya) = 1.8 x 10 -14 
erg s -1 cm -2 ), 2) within knot NW, < 50 pc (emission Lyaj v , 
F(Lya) = 5 x 10~ 15 erg s _1 cm -2 ), and 3) Northwest from knot 
SE and ~100 pc long (emission Lya v , F(Lya) = 1.6 x 10~ 14 erg 
s _1 cm -2 ). While Lya;;; shows a diffuse morphology, the Lya; v 
and Lya v contributions seem rather compact. However, the dif- 
fuse emission Lya;;; although dimmer than the other two com- 
ponents, is the strongest one when integrated over its large exten- 
sion , since it is observe d up to the northwestern extreme of the 
slit. lHaves et akl (120071) found also in Haro 1 1 that the diffuse 
emission surrounding the stellar clusters is the most prominent 
contribution to L(Lya). On the other hand, emission Lya; v is the 
weakest component, and unlike emissions Lya v and Lya^ it is 
fully located within one stellar knot. 



3.4. X-rays 

As explained in Sect. 12.21 the spectrum of the region labeled as 
NUCLEUS in Fig. |4]was analyzed. This region includes hard 
sources XI and X2, with the former completely dominating 
the hard emission of NUCLEUS. Fitting of the spectrum was 
performed with XSPEC going from simpler to more complex 
models, evaluating the statistical significance of the change at 
each step through the F-test and considering it valid if proba- 
bility significance was > 99%. Galactic absorption was fixed 
to Af(HI)c fl / = 6.3 x 10 19 crrT 2 , whereas the intrinsic absorp- 
tion was fixed to jV( HI),,, f = 7.0 x 10 19 cm~ 2 as derived by 
iLequeux et alj d 19951) . When freeing the hydrogen column den- 
sity of the intrinsic absorption it was found that the result of 
the fit and the value of N(HY)j nt did not change significantly, 
therefore we decided to fix it. Fo r the metallicity we assumed 
the value Z = Z Q /2 obtained by Davidgei dl989l) . The model 
which yielded a better fit is a composite of a hot plasma gas 
(hereinafter, HP, corresponding to the mekal m odel in XSPEC 
dMewe. Gronenschild. van den Oord et al.l ll985)) and a power- 
law emission (PL, in what follows) attenuated by both intrinsic 
and Galactic absorptions. HP would correspond to the diffuse 
soft emission filling the starbursting region and its surroundings, 
while the PL component is mostly dominated by XI. For the 
best fit achieved the temperature of the HP is kT = 0.71^q|q 
keV and the PL index F = l.S*!! - ?. No acceptable fit is obtained 
when considering only one single emitting component. Results 
assuming a Raymond-Sm ith model (raymond model in XSPEC 
(iRavmond & Smith|[T998l) ) for the HP instead of mekal yielded 
similar results. 
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Model 


kT (keV) 


r 


Metallicity 


V 


x 2 /y 




Norm" 


Norm" 


z 






1: wabs*zwabs*mekal 


0.81 




0.5 


16 


3.486 




2.9 x 10~ 5 










2: wabs*zwabs*zpowerlw 




2.1 




16 


2.947 






1.6x 10~ 5 








3: wabs*zwabs*(mekal+zpowerlw) 


71+011 


1 S+0-3 
l.S_0.4 


0.5 


14 


1.083 




1.6!° j x 10- 5 


9 + _\ x 10- 6 









a Units of normalization as in XSPEC. Hot plasma: 10~ 14 /{4n[D A (l + z)] 2 } jn e n H dV. Power-law: photons s -1 cm" 2 keV~' (1 keV). 



Table 9. X-ray fluxes and luminosities. 



Model 


F 0.4-2.4 keV " 

(erg s" 1 cm -2 ) 


^2.0-10.0keV a 

(erg s -1 cirr 2 ) 


I b 

-Lo.4-2.4keV 

(erg s" 1 ) 


I b (j PL ijHP \b 
LaO- 10.0 ke V (%4-2.4 keV ' ^0.4-2.4 keV > 

(erg s _1 ) 


3: wabs*zwabs*(mekal+zpowerlw) 


4.6 + _[° x 10- 14 


3.1^ x 10- 14 


2A + _° 5 6 x 10 39 


1.5!°;* x 1039 L2 -a5 



a Values of fluxes have not been corrected for neutral absorption. 
b Values of luminosities have been corrected for neutral absorption. 



Other options were tested, like fitting the spectrum with 
a two-temperature hot plasma, but large uncertainties were 
obtained for the higher temperature. Also, metal abun- 
dances were freed after modeling the HP with vmekal 
dMewe. Gronenschild. van den Oord et al ] fl985h . obtaining best 
fits when O and Mg abundances were Z(O)/Z(O) ~2.5 and 
Z(Mg)/Z(Mg) & ~4. However these models led to values of the 
reduced chi-squared of ^ 2 /v~0.7, which indicates an overpa- 
rameterization of the fitting a nd were therefore rejected. On 
the other hand, as explained by Ott, Fabian & Brinks (2005), a- 
elements are injected into the gas surrounding the stellar clusters 
by Type II SNe, whereas iron is mainly yielded into the interstel- 
lar medium (ISM) by Type I SNe, causing a time delay between 
both enrichments due to the different evolution of both types of 
stars producing each class of supernova. This causes a superso- 
lar value of the ratio Z(ff)/Z(Fe) in the gas inside bubbles blown 
by starbursts which have be en enriched by sup ernovae yields. 
This effect was modeled by Silich et al. (2001), who predicted 
an oversolar abundance of O compared to Fe during more than 
40 Myr after the onset of the starburst. In order to study the over- 
abundanc e of q-elements inside the hot gas we fo llowed the pro- 
cedure by Mar tin. Kobulnickv & Heckmanl d2002l) : the HP emis- 
sion was modeled by vmekal, He was fixed to the nebular abun- 
dance, Mg, Ne, Si and Ca were linked to the O abundance which 
was freed (considered to be the c-group, Z(a)), and the rest of 
the metals abundances were fixed to Fe abundance which was 
also left free (Fe group, Z(Fe)). However, no statistically signif- 
icant improvement in the fit was obtain ed. We studied the inclu- 
sion of S and Ar into the g -group as iGrimes et al. I (12001 and 
Kobu lnickv & Martini ((2010) did, but very similar results were 
obtained. 

The final results of the X-ray spectral fitting are shown in 
Table [8] Table [9] contains the X-ray luminosities for the best 
fitting model (Model 3), as well as the ratio of the contribu- 
tion of the PL component over the HP one in the soft X-ray 




1^ 
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Fig. 8. X-ray spectrum of Haro 2 with the best fitting model 
(Model 3 in Table. [8]). Thermal and power-law components are 
shown in dashed lines. Fitting residuals are also shown in the 
bottom panel. 



luminosity Lo.4-2.4keV- Finally Fig.[8]shows the X-ray spectrum 
of Haro 2 together with Model 3. Table [8] lists the values of 
the temperature kT of the HP, the PL index F, the metallicity 
of the gas assumed, the number of degrees of freedom v and 
the reduced chi-squa red y 2 /v. Errors corr e spond to a confi- 
dence level of 90%. IStevens & St rickland (Il998l) obtained a 
lower temperature value kT = 0.36 keV when analyzing the 
ROSAT/PSPC observation of Haro 2. They assumed a very low 
metallicity Z~Z / 10 and obtained a higher intrinsic hydrogen 
column density of A^(HI) = 2.2 x 10 21 crcT 2 with very high 
uncertainties, which might explain the disagreement. On the 
other hand the value we obtained in our study for kT is similar 
to the values found in other star-forming galaxies. For instance, 
Summers et al .1 (120041) obtained a value of kT - 0.75 keV for 
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the hotter component in their two-temperature composite hot 
plas ma when fitting the soft X-ray emiss ion of NGC 5253, 
and iMartin. Kobulnickv & Heckman] ([2002) obtained kT = 0.7 
keV for the component of the disk in the dwarf starbursting 
galaxy NGC 1569. lHaves et alJ d2007l) and lGrimes etal] d2007t) 
obtained kT = 0.69 keV and kT = 68 keV for Haro 11 
respec tively, andlGrimes et all (120051) and|Kobulnickv & Martin 
( 120101) calculated kT = 0.60 keV and kT = 0.63 keV for He 
2-10, respectively. However, these values are relati vely high 
when compared to the study by Strickl and et alJ d2004l) . who ob- 
tained kT = 0.11 keV and kT = 0.37 keV for a two-temperature 
model fitted simultaneously to eight star bursts whose X-ray 
point-sources had been removed. Also, iGrimes etalJ (120051) 
only found two sources with kT > 0.5 keV among seven dwarf 
starbursts. The kT value we propose for Har o 2 agrees better 
with the values found by IGrimes et alJ d2005l) for their sample 
ofULIRGs. 



IStrickland & Heckmanl (120071) found another source of hard 
X-rays in the star-forming nucleus of M82. After removing the 
point sources of the region they identified a diffuse hard X- 
ray component showing spectral features clearly different to the 
composite spectrum of the point sources removed and which ex- 
tends along the plane of M82. They identified this component as 
the wind fluid driving the superwind of the galaxy, emitting as 
much as i2.o-8.okev = 4.4x 10 39 erg s _1 , and produced by the col- 
limation of the massive stellar winds and SN ejecta of the newly 
formed masive stars. This emission can be spectrally fitted in the 
hard X-rays range either with a Bremsstrahlung thermal model 
(kT~3 - 4 keV), or with a power law with F > 2.0 yielding sim- 
ilar results. Nevertheless, we conclude that a wind fluid can not 
be responsible for the observed hard X-ray emission in Haro 2 
since the latter is unresolved, and it basically appears associated 
to XI. We rather argue that XI is most likely an ULX similar to 
the ones identified in other star-forming galaxies given the simi- 
larities in the values of its luminosity and photon index. 



As previously said, hard emission in Haro 2 is dominated by 
the point-like source XI. Attending to lPersic & Rephaehl (r2002) 
main contributors to the hard X-ray energy range 2-15 keV 
in starburst galaxies are High-Mass X-ray Binaries (HMXB) ac- 
creting gas onto a compact object. The gre a t spati al resolution of 
Chandra enabled Stric kland & Heckmanl (|2007) to extract the 
composite spectrum of the detected X-ray point sources in the 
nuclear starbursting region of M82 (excluding M82 XI identi- 
fied as an UltraLuminous X-ray Source). After fitting the hard 
energy range using a power law they obtained r~ 1.0. Similar 
values were found by IShaw Greeening et alJ d2009) for 18 ob- 
jects in M31 (T = 0.8-1.2). However, the fitting of Haro 2 X-ray 
spectrum points towards a photon index value of F~1.8 for XI. 
Even more, the X-ray luminosity of XI (io.4-10.okeV = 2.7 x 10 39 
erg s -1 ) is too large for a HMXB, and rather classifies it as 
an UltraLuminous X-ray Source (ULX), a type of hard-X-ray- 
emitting object with i2. 0-i0.0keV > 10 39 erg s" 1 and F values 
in the range 1.6 - 2.2 (Co lbert et al J I2004I) . Several of these 
sources have been highly studied in nearby galaxies where they 
can be sp atially separa t ed fro m the surrounding componen ts by 
Chandra. iKaaret et ail (l200lh and iMatsumoto et all (l200lh re- 
ported the observation of the ULX M82 X-l, a compact, hard 
X-ray emitter within the local starbursting galaxy M82 and 
170 pc far from its dynamical center. They found that its to- 
tal X- ray luminosity is Lx > 10 40 erg s -1 and Matsumoto et alJ 
( 200 1 ) calculated that the probability of this source being a back- 
ground AGN is < 1%. When fitted with an absorbed power law 
IKaaret. Simet & Lang I d2006l) ob tained a photon index F = 1 .67 
for this ULX. iKong et al.1 d2007l) analyzed another hard, com- 
pact source in M82 labelled as X42.3+59 only ~5" (-90 pc) 
away from M82 X-l. They also identified it as an ULX and they 
reported T = 1.3-1.7 andL^ o -io.oke v = (7.8- 1 l)x 10 39 erg s _1 . 
Even more, Stric kland et alJ d200ll) analyzed the X-ray emis- 
sion of the starburst galaxy NGC 3628, finding a compact object 
~20" (~1 kpc) away from the nucleus. They obtained that the 
model which best fitted the X-ray emission of the object was an 
absorbed power law, finding T =1.8 and Lo.3-8.okeV = 1.1 x 10 40 
erg s . IGrimes et al.l ([2007) also suggested that the hard, in- 
tense X-ray emitter they found analyzing a Chandra image of 
Haro 1 1 (also a local Lya-emitting starburst galaxy as Haro 2), 
could be an ULX, with a photon index F = 1 .6 and a luminosity 
of Lo.3-8.okeV = 2.2 x 10 39 erg s . Indeed, the X-ray luminosity 
Lo.4-10.OkeV = 2.7 x 10 39 erg s _1 and the photon index F~1.8 pro- 
vided by the X-ray spectral fit of Haro 2 point towards XI being 
an ULX, rather than an ordinary HMXB. 



4. Discussion 

The analysis of the observational data we have discussed in pre- 
vious sections indicates that the nucleus of Haro 2 has experi- 
enced very recent star formation, being dominated by a popula- 
tion of massive stars with ages in the range 3.5-5 million years 
and a total mass of ~4 x 10 6 M . The interstellar medium seems 
to be heavily distorted, with a patchy distribution of dust clouds 
which blocks most of the UV continuum emitted by a high frac- 
tion of these stars, which are identified only by their contribution 
to the ionization of the gas and the heating of the dust parti- 
cles. Moreover, while Haro 2 is surrounded by a relatively dense 
cloud of neutral gas, parts of it have been accelerated by the re- 
lease of mechanical energy by these powerful starbursts, opening 
kinematical holes which allow the escape of Lya photons out of 
the ionized regions, albeit in the form of P Cyg profiles. The 
release of mechanical energy has also led to powerful, diffuse 
soft X-ray emission around the central star-forming knots. Even 
more, Haro 2 presents a strong hard X-ray emission of the order 
of the soft X-ray emission, that is dominated by a single, unre- 
solved source which is an ULX candidate. Finally, an extended 
and diffuse component of Lya emission has been detected whose 
origin is unclear. This diffuse emission is spatially not correlated 
either with the UV continuum, nor even with the Balmer lines 
emission or Balmer decrement, but on the other hand it might 
originate close to the soft X-ray emitting regions. 

4. 1 . Past star formation episodes 

Previous episodes of massive star formation should have con- 
tributed to the chaotic structure of gas and dust in the nucleus of 
Haro 2. Some authors claim indeed that star formation has taken 
place in Haro 2 previously t o the current starbursting episodes. 
Fanell i. O'Connell & Thu an (1988) studied signature lines of 
stars of different spectral types in spectra taken by the IUE. They 
found that the presence of 07-BO V stars in Haro 2 reveals a 
young burst in the galaxy having ocurred within the past 10 Myr. 
However, due to the presence of A-type stars, they also argued 
that Haro 2 must have experienced two previous starburts, the 
most recent having ended less than 20 Myr ago. From this data, 
Summers, Stevens & Strick land! d200lh claimed that the star for- 
mation history of Haro 2 is composed of a recent burst with age 
~5 Myr, another of 20 Myr, and a previous one which took place 
around 500 Myr ago. Also, they found that the older the burst, 
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the larger the number of stars produced, and hence they con- 
cluded that the intensity of the episodes of the star formation 
must have decreased with time, maybe due to the lack of gas 
available for star production. 

In order to check t he presence of past sta r-forming episodes 
we used the models bv lMartfn-Manion et alJ (120081) which com- 
bine codes of evolutionary population synthesis models, chem- 
ical evolution and photoionization to model the intensity of op- 
tical nebular lines, magnitudes and abundances of metals in a 
star-forming region. The models assume a history of star for- 
mation composed of bursts equally-separated, attenuated in time 
and corrected for the efficiency of the stellar production from the 
available gas. The production of stars from gas which has been 
chemically enriched by stellar winds and supernovae is consis- 
tently considered, as well as the effect of this enrichment on the 
intensity of nebular lines produced by photoionization. Our aim 
was to use the models to check whether previous star forma- 
tion episodes are needed to reproduce t he integrated intensity of 
the em ission nebular lines measured by Mo ustakas & Kennicuttl 
(2006). Lines obs erved in the spectrum and whose intensities 
are predicted by iMartfn-Manion et alJ d2008l) are: 0[II]3727, 
O[III]5007, O[I]6300, Ha, N[II]6584, S[H]6716 and S[II]6731. 
Intensity is scaled to H/3 and observed lines were corrected for 
stellar absorption and Galactic and intrinsic extinction. Since 
L(Ha)/L(Hj8)~2.86 was assumed to correct for internal extinc- 
tion we did not include Ha in the analysis to avoid biases to- 
wards models showing similar values. 

Given the integrated value measured by 
iMoustakas & Kennicuttl d2006l) for Haro 2 £W(H/3)=44 A, 
we only considered models with predicted EW(H/3) in the 
range 34 A<EW(H/3)< 54 A. Models consider bursts equally- 
separated in time by either 0.05, 0.1, 0.5 and 1.3 Gyr. Results 
show that models which best reproduce the intensity values 
of the lines are those 1) having had previous star formation 
bursts, either each 0.5 Gyr or 1.3 Gyr, and 2) last burst 
having taken place 3-5 Myr ago. Point 2) agrees with 
the ages of the bursts w e have calculated, whereas point 1) 
reinf orces the idea by iFanelli. O'Connell & Thuanl dl988h 
and ISummers. Stevens & Strickland! (120*01 ) that intense star 
formation in Haro 2 took place as much as hundreds of millions 
of years before the current episode. Therefore these results 
confirm that previous star formation in Haro 2 occurred and 
is responsible for some of the features observed in the optical 
spectrum. 



4.2. X-rays 

Soft X-ray emission in a starburst is expected to be produced by 
1) the heating of the diffuse gas surrounding the stellar clusters 
up to temperatures of millions of Kelvin, due to the stellar winds 
and supernovae injecting mechanical energy into the medium, 
and 2) by the emission of supernova remnants (SNRs) during the 
adiabatic phase. The contribution by the direct emission of indi- 
vidual starsjjsjiejyigjbje_wjie^ the previous compo- 
nents (Cervino, Ma s-Hesse & Kun th 2002). Evolutionary popu- 
lation synthesis models CMHK02 predict the soft X-ray emis- 
sion of a starburst after defining the physical properties of the 
burst. CMHK02 model the component 1) through a Raymond- 
Smith plasma with kT = 0.5 keV, considering that a fraction e xe s 
from the mechanical energy is finally converted into X-ray lumi- 
nosity. Typ ically, values e XK ff= 1-10% are found for star-forming 
objec t s (ISummers. Stevens & Stric kland 2 0011: I Summers et alJ 
2004; iMas-Hesse. Otf-Floranes & Cervino! 120081) . On the other 



hand, the soft X-ray emission by SNRs is considered to be 
due to the reversed shock which heats the gas inside the shell 
and which is modeled by a composite of three Raymond- 
Smith plasmas with kT = 0.23, 0.76, and 1.29 keV. See 
ICervino. Mas-Hesse & Kunthl d2002l) for a more detailed de- 
scription. 

As already explained, using the ROSAT/PSPC observation of 
Haro 2 lStevens& Strickland! (119981) obtained values kT = 0.36 
keV for the temperature of the hot plasma and A^(HI) = 2.2x 10 21 
cm" 2 for the absorption, which yielded a higher soft lumi- 
nosity Lrqsat~1-4 x 10 40 erg s" 1 wh e n the y assumed Z ~ 
Z o /10. ISummers. Stevens & Strickland! (l200ll) argued that the 
high u ncertainties obtained for Af(HI) by IStevens & Strickland! 
(1998) caused also large uncertainties in the X-ray luminos- 
ity, which was constrained in the range Lrosat~2 x 10 39 - 
1.4 x 10 42 erg s" 1 and took L R ns AT ~2 x 10 39 erg s _1 as valid. 
Sum mers. Stevens & Strickland (2001) estimated the injection 
of mechanical energy into the medium per unit of time Ek by 
two meth ods: 1) applying the mode l for s tarbur st-driven out- 
flows by ICastor. McCrav & Weaver! (1 19751) and IWeaver et alJ 
(1 19771) assuming that the radius of the superbubble blown by 
the star-forming activity in Haro 2 is Rs = 1.1 kpc, and 2) using 
the ratio Ek/Li, i from SB99 models, where Lb \ was considered 
to be well represented by the observed far infrared luminosity 
£ftr= 1.4 x 10 43 erg s _1 . After performing these calculations for 



^FIR = 

several ages, they assumed the value Ek- 1 .2 x 10 41 erg s _1 con- 
sidering an age of the burst of 5.8 Myr. Thus, with their values of 
Ek and LpncaT- ISummers. Stevens & Strickl and (2001) obtained 
e xe ff~0.02 for Haro 2. 

We calculated the X-ray luminosity predicted by the 
CMHK02 models for a composite of two starbursts with the ages 
and masses found for Population 1 and Population 2 in Sect. 
I3.2l and included in Table [7] Metallicity of the models was fixed 
to Z = 0.008 since it is the closest value to the observed one 
Z~Z /2 (lDavidgeHl989l) . As shown in Fig. [8] soft X-ray emis- 
sion is not strictly due to gas heated by the starburst activity 
(HP component), but also by the XI emission (PL component). 
Actually, as indicated in Table|9]soft X-ray emission due to HP is 
~45% of the total soft X-ray emission observed, i.e. gas heated 
by the starburst activity emitts Lo.4_2.4keV~l-l X 10 39 erg s . 
CMHK02 can reproduce this soft X-ray luminosity value for 
the composite of the bursts explained assuming e xe jf= 0.004, 
i.e. the efficiency in the conversion of mechanical energy in- 
jected into the ISM by SNe a nd stellar winds in Haro 2 into 
soft X -ray luminosity is ~0.4%. Summ ers. Stevens & Strickland 
(2001) found the much higher value of 2% since they did not 
subtract the emission of XI from the total observed luminosity. 
If we consider that the whole total luminosity Lq .4-2.4 keV ~2.4 x 
10 39 erg s~' is due to gas hea t ed by the starburst activity as 
ISummers. Stevens & Strickland! (l200lh did, we obtain an effi- 
ciency of ~1% which is in better agreement with their value. 
When calculating the efficiency e xe jf for galaxies where no sub- 
traction of the emission from compact, point-like emitters is pos- 
sible (due to lack of spatial resolution, for instance), then com- 
parison must be made with efficiency values obtained with the 
total soft X-ray lumino sity. As explained, this is the c ase fo r 
the value reported by Summers. Stevens & Strickland! d200ll) . 
Given the uncertainties in the calculation of the X-ray lumi- 
nosity predicted by the CMHK02 models and in the estimation 
of the contribution by XI to the soft X-ray emission, the range 
e xe ff= 0.004 - 0.01 should be assumed for the value of this effi- 
ciency in Haro 2. 



13 



Otf-Floranes et al.: Multiwavelength analysis of starburst in Haro 2 



4.3. Lya emission 

Starbursts in the nucleus of Haro 2 dominate the ultraviolet 
emission of the galaxy (mainly by massive stars in dust-free 
environments), as well as its ionizing power (by the most red- 
dened young stars), which produces conspicuous nebular emis- 
sion lines. The spatial profiles of the UV continuum and the Lya 
and Ha emission lines along directions Dl and D2 are shown 
in Fig. [7] together with the observed ratios Lya /Ha and Ha/H/3. 
Also, the profile of the expected Lya /Ha ratio is included, which 
was calculated assuming case B recombination with T e = 10 4 
K and n(e) = 500 cm" 3 (Lya/Ha = 8.7, Ha/H/3 = 2.87), ap- 
plying the expected internal dust reddening as derived from the 
observed Ha/H/3 ratio, as well as the corresponding Galactic ex- 
tinction. Values of signal-to-noise ratio larger than 3 were ob- 
tained in the profiles by lowering the spatial resolution of the 
emission lines data when needed. Specifically, the resolution 
in Ha/H/3 was lowered down to 2.5 arcsec pixel 1 , whereas 
Lya/Ha was calculated for 0.5 arcsec pixeL 1 and 1.0 arcsec 
pixeL 1 for directions Dl and D2, respectively. As seen in the 
figure, the profiles of the UV continuum and the Ha emission 
line are tightly spatially correlated, as expected since the B aimer 
lines are mainly produced by the recombination of gas ionized 
by the young massive stars. 

On the other hand, Lya is not always spatially coupled to 
the stellar continuum, Balmer lines emission or Ha/H/3 ratio. 
While compact Lya components Lya;, Lya; v and Lya v are lo- 
cated within (or close to) a Balmer- line-emitting region, diffuse 
emissions Lya;; and Lyani appear in regions where the contin- 
uum and other lines are very weak. In what follows we will an- 
alyze these 2 sets of regions independently, since their origin 
is probably different. Lyai, Lyaj v and Lya v are apparently pro- 
duced by recombination in the gas ionized by the massive stellar 
clusters, but the variable conditions of the neutral gas surround- 
ing these regions (mainly kinematics and column density) deter- 
mine whether the emission can escape or not, and whether the 
line profile is kinematically distorted. The final result is the al- 
ternance of peaks and valleys observed in the Lya spatial profile 
between knots SE and NW along slit D2 (see Fig. |7). We want 
to stress that this structure is not reflected in the spatial profile of 
Ha since this line is not affected by resonant scattering. 

The Lya emission identified in regions Lyau and Lya;;; 
shows strikingly different properties. First of all, the Lya/Ha 
ratio is much higher in these areas than close to the stellar clus- 
ters (see Fig. |7]l. Ratios as high as Lya/Ha = 3 and Lya/Ha = 6 
are reached along directions Dl and D2, respectively. These ob- 
served ratios are much larger than predicted applying case B and 
the reddening derived from the Balmer decrement, representing 
Lya emission ~3 and >100 times higher than expected, respec- 
tively. Furthermore, it is remarkable that along D2, Lya/Ha is 
highest where Ha/H/3 also reaches its highest value (NW end of 
the slit), leading initially to reject the hypothesis that the rela- 
tively large Lya/Ha values in the diffuse regions are originated 
simply by a low abundance of dust. 

Diffuse Lya emission similar to Lya;i and Lya^ has been 
found in other local Lya-emitting galaxies. lHaves et"a"fl (2005) 
found that 70% of the total Lya radiation of the galaxy ES0338- 
IG04 is diffuse emission located in a halo surrounding the cen- 
tral regions and lac king strong stellar continuum. Furthermore, 
lHaves et"al] J2007) concluded that 90% of the Lya emission in 
Haro 1 1 is also diffuse and is located in a featureless halo, where 
Lya flux is much higher than expected fro m recombination , 
when scaled from the observed Ha intensity. lAtek et aTl ([2008 ) 
studied local values of the Lya/Ha ratio in a sample of local 



Ly a-emitting s t arburs t gal axies. Together wi th the values found 
bv lHaves etafl (120051) and lHaves et aD (120071) for ESO338-IG04 
and Haro 1 1 respectively, they rep orted those for IRA S 0833 and 
NGC 6090. For the latter sources lAtek et al J d2008l) also found 
Lya emission values higher than predicted from case B recom- 
bination considering the Ha fluxes observed, when applying the 
extinction derived from the Balmer decrement (Ha/H/3). 

These properties indicate that the origin of the Lya emis- 
sion in the diffuse regions is apparently different from that of the 
compact knots dominated by the massive clusters. We postulate 
that the Lya emission in these diffuse regions could be enhanced 
with respect to Ha either by an additional emission component, 
associated to the hot plasma identified in the X-ray image, by 
leaking of the Lya photons created in the ionized regions, af- 
ter multiple scattering within the neutral clouds, by reflection 
on a dust screen, or by a combination of the 3 processes. In the 
next subsections we describe in some detail these different, non- 
exclusive models. 

4.3.1. Photoionization by the hot plasma 

For this model we assumed that the Lya diffuse emissions Lya;; 
and Lya;;; in Haro 2 have been produced by ionization domi- 
nated mainly by the hot plasma responsible for the intense soft 
X-ray radiation, and not by the massive stellar clusters as for 
Lyai, Lyaiv and Lya v . As we observe in Fig. [4] there seems to 
be some spatial correlation between the diffuse Lya components 
and the diffuse, soft X-ray emission, originated by the heating in 
shock fronts of the interstellar gas surrounding the stellar clus- 
ters. We postulate that this X-ray radiation might have ionized 
the nebular gas in these regions, in an environment far from case 
B recombination conditions. 

To check this possibility we used t he photionization code 
CLOUDY 10.00 dFerland et all 119981) . Ratios Ha/H/3 and 
Lya/Ha were calculated in two extreme cases in which a spheri- 
cal cloud of gas is ionized by a central source, which is a massive 
star in Model 1, and a hot plasma in Model 2. In the former we 
modeled the continuum with a Kurucz star with T~4- x 10 4 K 
and emitting Q(H) = 10 50 s _I ionizing photons, whereas the 
electronic density was assumed to be «(e) = 10 3 cirT 3 . Given a 
typical distance between the star and the cloud of gas of ~1 pc, 
the value of the ionization parameter is U = 0.1. On the other 
hand, for Model 2 we assumed a Bremsstrahlung emission with 
T = 10 7 K and n(e) = 10 3 cm" 3 , keeping the ionization pa- 
rameter within U = 10~ 8 - 10~ 5 . This extreme conditions were 
chosen to reproduce the physical properties of high-density fila- 
ments compressed by the shock fronts. 

The emission line ratios computed for Model 1 are 
(Ha/HyS) Mo dei 1 = 2.9 and (Lya/Ha)Modeii = 10.1, indeed very 
similar to the predictions of the standard case B, as expected. On 
the other hand, for Model 2 we obtained (Ha/H/T)Modei2 in the 
range 3.0-4.0, and (Lya/Ha)Modei2 from 7.0 to 12.0. It is impor- 
tant to note that Model 2 predicts systematically higher Ha/H/3 
values than case B. This implies that if the contribution by the 
hot plasma is important, we could overestimate the abundance 
of dust if derived assuming case B conditions. 

We would expect a physical scenario with a smooth tran- 
sition from Model 1 to Model 2 conditions as we move away 
from the vicinity of the massive clusters to more external re- 
gions. Fig. [7] shows compact Lya emission, if any, close to the 
stellar clusters in Haro 2. Assuming that Model 1 conditions pre- 
vail in these regions, reddening correction as derived from the 
case B Ha/H/3 ratio should be applied. The Lya/Ha ratios ob- 
served around these regions are generally lower than expected 
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from case B, even taking into account the effect of reddening. 
This is especially noticeable around Lya n and Lya v . We know 
already from UV spectroscopy that scattering in the outflow- 
ing neutral gas surrounding knot SE, together with the relatively 
high dust abundance (as derived from the observed Ha/H/3 ra- 
tio), are responsible for the low escape fraction of the Lya pho- 
tons originated on this area. 

On the contrary, the diffuse Lya emission is originated in re- 
gions spatially separated from the stellar clusters, where the ion- 
ization might be dominated by the hot plasma. Since Model 2 
predicts higher intrinsic Ha/Hfi ratios than Model 1, the extinc- 
tion derived in the diffuse emitting regions would be much lower 
than assuming case B, allowing for higher Lya escape fractions 
in these areas. The ratio measured on Lyqjj, Lyq/Hq< 3.0, i s 
compatible with the results found by Mas -Hesse et alJ (|2003). 
They found that the diffuse Lya emission exhibits a clear P Cyg 
profile all along slit Dl, at any distance from the stellar knot 
SE. This implies that a significant fraction of the emitted Lya 
photons (the blue wing of the line) has been affected by reso- 
nant scattering, decreasing the outcoming line flux in this area. 
On region Lyajjj, both the observed Ha/H/3 and Lya/Ha ratios 
could be consistent with the presence of shock-heated filaments, 
with small or none attenuation by dust. Indeed, the strong X-ray 
flux in this area would most likely destroy most of the dust parti- 
cles. The fact that the Lya/Ha ratio is very close to the intrinsic 
predictions by Model 2 requires also that the medium in front 
of this region is transparent to Lya photons: either the column 
density of neutral gas is small enough, or it is outflowing with 
a projected velocity large enough not to affect the Lya intensity 
(v ou t S 400 km s _1 ). Unfortunately the spectral resolution is not 
enough to analyze the spectral profile of Lya^ . 

Model 2 predicts L(Ha)/Lo.4-2.4keV~ 10~ 3 - 10~ 2 , where 
L(Ha) and Lo.4-2.4keV are the intrinsic Ha luminosity produced 
and the ionizing soft X-ray luminosity emitted by the hot 
plasma, respectively. We measured L(Ha)/Lo.4-2.4keV in both 
the regions dominated by the stellar clusters and those show- 
ing diffuse Lya emission. For the compact regions we found 
L(Ha)/Lo.4-2.4keV~300, while in the diffuse regions, far away 
from the stellar clusters, L(Ha)/Lo.4-2.4keV~l- This severe de- 
crease in the L(Ha)/Lo.4-2.4keV ratio when moving towards re- 
gions showing diffuse emission reinforces the idea of ionization 
in these sites being dominated by the hot plasma itself. Values 
of L(Ha)/Lo.4_2.4keV~10 - 100 are measured in the intermediate 
regions, where the ionizing flux of the massive stars increases 
with respect to the hot plasma one as we move closer to the stel- 
lar clusters. Although the measured L(Ha)/Lo.4-2.4keV is always 
higher than the absolute values predicted by Model 2, we should 
be aware that the hot plasma filaments which are expected to ion- 
ize the local surrounding medium can not be resolved, and hence 
the flux values and ratios obtained are average values. Measuring 
L(Ha)/Lo.4-2.4keV locally would require a higher spatial resolu- 
tion and better statistics in the X-ray observation. However, the 
observed trend in L(Ha)/Lo.4-2.4keV does point towards the dif- 
fuse Lya emission being produced by the ionization caused by 
the hot plasma, whereas the compact Lya emission is consistent 
with being originated in the gas ionized by the massive stars. 

Finally, while the simulations with CLOUDY suggest that 
the hot plasma could dominate the Lya emission in the exter- 
nal regions of Haro 2, we have not been able to reproduce all 
the spectral features observed in the regions showing Lya;; and 
Lya;;; emissions. Whereas the ionizing hot plasma models, com- 
bined with some attenuation, can reproduce the observed ra- 
tios Lya/Ha, Ha/H/? and L(Ha)/Lo.4-2.4keV, the predictions for 
some other lines (like Fe II 4300 and Fe II 6200) are overesti- 



mated when compared with the observations. Furthermore, other 
lines present in the spectra such as 0[II]3727 and 0[III]4959, 
5007 are underestimated by these models. Nevertheless, we want 
to stress that our models are based on a simplistic scenario 
assuming spherical geometry, while X-ray-emitting filaments 
should be significantly more complex in shape, structure and 
physical conditions. A more detailed modelling of the geometry 
of these filaments is out of the scope of this work. We want to 
retain here only the idea that a very hot, X-ray-emitting plasma 
could generate an additional Lya emission component under 
conditions that would differ significantly from case B ones, and 
which might be the major contributor to the diffuse Lya emission 
in Haro 2. 



4.3.2. Scattering by the neutral gas 



As we have already discussed, Lya photons are 
resonantly scattered by neutral hydrogen atoms. 
IVerhamme. Schaerer & Masellil d2006l) have modelled the 
properties of Lya radiation transfer through a neutral medium, 
assuming different conditions of dust abundance, kinematics of 
the gas, column density,... In the absence of dust, Lya photons 
could travel long distances through the neutral cloud, before 
leaking from its surface, escaping from the galaxy. Nevertheless, 
the multiple scattering would at least modify their distribution 
in energy, leading to a two-peak emission profile with almost no 
emission at the central wavelength. If the gas has some velocity 
with respect to the ionized region where the Lya photons were 
produced, a P Cyg-like profile would be expected. The multiple 
scattering increases the mean free path of the photons through 
the neutral cloud, increasing the probability of being absorbed 
by dust. 

The presence of neutral gas around the Haro 2 nuclear star- 
burst is evident from the observations. The P Cyg profile of 
the Lya emission line along slit Dl indicates that a least knot 
SE is surrounded by outflowing neutral gas. The regions where 
the Lya emission disappears completely along slit D2 are most 
likely also covered by neutral gas clouds, in different kinemati- 
cal states. Therefore, it might be possible that the Lya photons, 
formed in the ionized regions close to the stellar clusters, but be- 
ing scattered by these neutral clouds, end up leaking at relatively 
large distances. Since Ha would not be affected by this process, 
the net result would be an enhancement of the Lya emission on 
extended areas far away from the ionized regions, as proposed by 
lHaves et"a"D (120071) . which in principle would be consistent with 
the observations. However, in that case the Balmer lines would 
be produced by the ionizing flux originated in the stellar clus- 
ters, and conditions similar to case B would apply. The observed 
Ha/Hy6 values would indicate a high abundance of dust particles 
in front (or mixed with) the clouds surrounding Lya;; and Lya^ 
regions, and they would most likely destroy most of the scat- 
tered Lya photons. Nevertheless, while the Lya photons would 
be leaking from the surface of the neutral clouds, the residual 
Balmer lines emission in these regions would probably be origi- 
nated at a different location, so that the relatively high extinction 
derived from the Ha/H/3 ratio would not apply to the Lya emis- 
sion. A detailed analysis of the Lya spectral profile of the Lya^ 
region, not possible with the low-resolution existing data, would 
be required to compare it with the predictions of radiation trans- 
fer models. 
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4.3.3. Dust reflection 

If we do not consider the possible contribution to ionization by 
the hot plasma, and assume that conditions close to case B apply 
over most of the nuclear region of Haro 2, then the high Ha/Hfi 
ratios (especially on region Lya;;;) would point to the presence 
of dense dust clouds or filaments at certain distance from the 
stellar clusters. These dust clouds could be reflecting the emis- 
sion originated close to the starburst knots. Since UV photons 
are more efficiently dispersed than optical photons, an enhance- 
ment of the Lya emission over the Balmer lines in these regions 
would be naturally produced, explaining the higher Lya/Ha ra- 
tios in these extended regions. In principle we would expect also 
Ha/H/3< 2.86 in these regions, but since dispersion would be so 
much efficient for Lya wavelengths than for Ha or HB, it might 
be well that the observed Ha and HB emissions are dominated 
by the component originated locally, and therefore attenuated by 
the local dust. The local Lya would be almost totally destroyed, 
as expected from the high E(B-V) derived. But the reflected 
Lya component (scattered at the surface of the dust clouds, and 
therefore not affected by the internal extinction) could be strong 
enough to produce the observed Lya/Ha values. The main flaw 
of this model is that it predicts lower Lya equivalent widths than 
observed, since the UV continuum would also be reflected. 

These different mechanisms are not mutually exclusive, and 
each of them could be contributing to a different extent to 
the Lya emission d etected in starburst galaxies like Haro 2. 
IScarlata et alJ (120091) proposed indeed a scenario which com- 
bines both the scattering by neutral gas and the reflection by 
dust clouds, based on a clumpy distribution of neutral clouds 
rich in dust, surrounding the ionized region. While the Balmer 
lines and the UV continuum would be attenuated by the dust 
particles along the line of sight, scattering and/or reflection at 
the surface of the neutral clouds would allow the Lya photons to 
find their way through the medium with a rather high efficiency, 
leading to high Lya/Ha ratios when observed globally. They 
argued that this scenario could explain the integrated values of 
Lya/Ha =1-5 measured in a sample of GALEX Lya-emitters 
at z = 0.18 - 0.28. Mapping the distribution of the dust clouds 
in Haro 2 with very high spatial resolution would be needed to 
verify if this scenario could explain the diffuse Lya emission. 

However, we argue that ionization by the hot plasma is 
the dominating source responsible for the diffuse emission in 
Haro 2. Despite the limitations of this model concerning the pre- 
dicted intensities of some emission lines, it can naturally explain 
the spatial coupling of the Lya;; and Lya;;; components with the 
soft X-ray emission, the measured values of the Lya/Ha and 
L(Ha)/Lo.4-2.4keV ratios and the high escape fraction of Lya pho- 
tons from regions with high Ha /H/3 ratios. 

If the diffuse Lya emission detected in Haro 2 and other star- 
burst galaxies is originated by ionization by hot plasma, and not 
directly by the massive stars, the use of Lya could overestimate 
the star formation rate in galaxies at hight redshift b y a signifi- 
cant factor (w hen integrated over the whole galaxy. lHaves et a0 
( 20051 120071) found that the diffuse emission is a dominant frac- 
tion of the total Lya emission observed). This adds to the diffi- 
culty of determining the escape fraction of Lya photons orig- 
inated dir ectly by stellar i onization along the history of the 
Universe (lHaves et al.ll201 lh . A further correction should be ap- 
plied to the procedure used to convert Lya luminosity into num- 
ber of ionizing photons, and thus star formation rate. Calibrating 
the strength of these corrections requires a detailed analysis of 



a larger sample of Lya-emitting starburst galaxies in the local 
Universe. 

5. Conclusions 

We have carried out a complete multiwavelength spectral and 
photometric analysis of the star-forming nucleus of the local 
Lya-emitter Haro 2, focusing on the Lya emission and its rela- 
tion with different starburst-related parameters. UV, optical, NIR 
and X-ray images, as well as UV-optical spectral observations 
with high spatial and spectral resolution, have been analyzed. 

- We have characterized the starburst of Haro 2 by compar- 
ison with evolutionary population synthesis models, repro- 
ducing the UV, Ha and FIR luminosities values as well as 
the properties of the SilV and CIV stellar absorption lines. 
We have identified a population of young massive stars (ages 
-3.5 - 5.0 Myr) with a total stellar mass M~4 x 10 6 M (for 
an IMF with mass limits of 2 - 120 M ), and whose con- 
tinuum is affected by a differential extinction with E(B-V) 
values in the range E(B-V)=0.035- > 0.5. Stars located in 
dust-free regions dominate the UV emission, whereas those 
within dust-rich clouds are completely obscured in the UV, 
but are the main contributors to the Ha and FIR luminosities 
of the galaxy. 

- We have identified a diffuse soft X-ray emission extended 
over the whole nucleus. This component is attributed to the 
gas heated by the release of mechanical energy by the stel- 
lar winds and supernovae from the starburst, as well as to 
supernova remnants. We have estimated that ~QA - 1% of 
the mechanical energy injected is converted into soft X-ray 
emission, which is a typical value for this kind of sources. 
A hard, unresolved X-ray emission is located on the SE star- 
forming knot, apparently originated by an UltraLuminous X- 
ray (ULX) source related to the starburst episode. 

- Both compact and diffuse Lya emission components, to- 
gether with regions of total absorption, are observed along 
the major and minor axes of the star-forming nucleus of 
Haro 2. Intensities of the compact Lya emission are lower 
than those expected from Ha fluxes and case B recombina- 
tion theory, considering dust extinction. On the other hand, 
the diffuse Lya emission is stronger than expected from Ha, 
and extends over 850 pc and > 600 pc along minor and ma- 
jor axes, respectively. We have found that the Lya emission 
is spatially decoupled from the UV continuum, the Balmer 
lines emission and the Ha/H/? ratio, but the diffuse Lya com- 
ponent is apparently correlated with the diffuse soft X-ray 
emission. 

- Outflowing of the neutral hydrogen surrounding the stellar 
clusters, energized by the starburst activity, enables Lya pho- 
tons from both compact and diffuse components to escape 
from Haro 2, leaving the signature of a clear P Cyg profile 
in the emission line. Nevertheless, the irregular distribution 
and kinematical properties of the expanding neutral gas lead 
the Lya photons produced in recombination sites to escape 
only in certain regions. Ha photons are not affected at all by 
the kinematics of the neutral gas, and escape directly from 
all ionized regions. This explains the spatial decoupling be- 
tween the Lya and Ha compact emissions. 

- Compact Lya emission is produced by recombination in the 
gas ionized by the massive stellar clusters. The combined 
effect of resonant scattering by neutral gas and attenuation 
by dust, as traced by the P Cyg profiles and the high values 
of Ha/Hy6 ratio respectively, decreases the escape fraction of 
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the Lya photons and yields observed fluxes below the case 
B recombination predictions. 

- On the other hand, the diffuse Lya emission has apparently a 
rather different origin. We argue that the shock fronts heated 
to X-ray emitting temperatures by the release of mechanical 
energy could ionize the surrounding nebular gas under con- 
ditions far from case B recombination, as suggested by the 
spatial correlation between the soft X-ray emission and the 
diffuse Lya component. Simple modelling with CLOUDY Re- 
produces indeed the Lya/Ha, Ha/H/3 and L(Ha)/Lo.4-2.4keV 
ratios measured in this region. 

- Nevertheless, our results do not reject other physical pro- 
cesses proposed previously that could also contribute to en- 
hance the intensity of the diffuse Lya emission line. Lyar 
photons, though originated close to the massive stars, could 
be leaking from the surface of distant neutral clouds after 
multiple scattering within the gas, depending on its distribu- 
tion and the abundance of dust. Moreover, under certain cir- 
cumstances the dust clouds could even contribute to enhance 
the diffuse Lya emission by differential reflection. 
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Fig. 7. Left: slit along direction Dl. Top-left: HST/STIS G140M spectral image of Haro 2. Middle-left: Emission profiles of Lya 
(blue) and UV continuum (1 175 A rest frame, black) from STIS data, and Ha (red) from WFPC2 data. Bottom-left: Spatial profiles 
of Lya (blue) and observed Lya/Ha (black) and Ha/H/? (from WHTflSIS data, in red). Right: slit along direction D2. Top-right: 
HST/STIS G140L spectral image. Middle-right: Emission profiles of Lya (blue) and UV continuum (1227 A, black) from STIS 
data, and Ha (red) from WFPC2 data. Bottom-right: Spatial profiles of Lya (blue) and observed Lya/Ha (black) and Ha/H/3 (from 
WHTflSIS data, in red). The geocoronal Lya line has been blackened in both spectral images at the top. In the middle and bottom 
panels positive x-axis corresponds to Southwest (left) and Southeast (right). Vertical scale in the middle panels corresponds to the 
Lya and Ha fluxes. To ease the comparison, the dashed line shows the Ha profile scaled by a factor 0.4 (left) and 0.085 (right). The 
vertical scale in the bottom panels corresponds to the Lya/Ha and Ha/H/? ratios. The brown dashed lines represent the expected 
Lya/Ha ratio, assuming case B recombination and applying the internal reddening derived from the observed Ha/H/3 ratio (Balmer 
decrement), as well as the Galactic extinction. 
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